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Abstract 

We present Monte Carlo Simulations using an equilibrium lattice- 
gas model for the electrosorption of CI on Ag(100) single-crystal sur- 
faces. Fitting the simulated isotherms to chronocoulometry experi- 
ments, we extract parameters such as the electrosorption valency 7 
and the next-nearest-neighbor lateral interaction energy nnn . Both 
coverage-dependent and coverage independent 7 were previously stud- 
ied, assuming a constant 4> nan [I. Abou Hamad, Th. Wandlowski, 
G. Brown, P.A. Rikvold, J. Electroanal. Chem. 554-555 (2003) 211]. 
Here, a self-consistent, entirely electrostatic picture of the lateral in- 
teractions with a coverage-dependent 4> nan is developed, and a rela- 
tionship between nnn and 7 is investigated for CI on Ag(100). 

Keywords: Chlorine electrosorption; Lateral interactions; Electrosorp- 
tion valency; Chronocoulometry; Continuous phase transition; Lattice-gas 
model; Monte Carlo simulation. 
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1 Introduction 



Studies of lateral interactions between adsorbed particles are motivated by 
the need to understand the origin of the wide variety of ordered overlayers and 
phase transitions at fractional adsorbate coverage on metal surfaces. These 
interactions have contributions ranging from short-range and van der Waals 
to long-range dipole-dipole, and lattice-mediated interactions jl . Hard- 
square short-range interactions and dipole-dipole long-range interactions are 
the major contributions to the lateral interactions for bromine adsorption on 
Ag(100) 2J. In this paper, we explore the validity and applicability of such 
a model for adsorption of chlorine on Ag(100). 

Halide electrosorption on single- crystal metal electrode surfaces is a good 
model system for studying the properties of the electrode-electrolyte interface 
in an electrochemical cell. Due to its relative simplicity, it can be used to dis- 
tinguish between the various contributions of different parameters according 
to the effect of their inclusion on the overall behavior of the system. A mean- 
field approach is not sufficient, even for the description of one of the simplest 
halide-electrosorption systems Br/Ag(100). However, a simple lattice-gas 
model with constant parameters is sufficient to describe its equilibrium [21 E] , 
and dynamic [3] properties. While the electrosorption of Br on single- crystal 
Ag(100) from aqueous solution has been extensively studied as an example 
of adlayer formation in an electrochemical system [21 El El El E] , less attention 
has been given to the electrosorption of CI [3 El ED] on Ag(100). A lattice- 
gas model with constant parameters is not sufficient to describe Cl/Ag(100), 
therefore this system can be used to further investigate the nature and char- 
acteristics of the lateral interactions between the adsorbed halide atoms. In 
particular, we here develop a self-consistent picture of variable lattice-gas 
parameters based on the resident charge on the adatoms being coverage de- 
pendent or electrode-potential dependent (through the coverage). 

The rest of this paper is organized as follows. In Section El we describe an 
electrostatic model of the adlayer that is used in the simulations, the lateral 
interaction energies, and the Monte Carlo methods used. A brief description 
of the experimental procedure is given in Section El The results of fitting the 
simulations to experimental data are detailed in Section EJ followed by a brief 
comparison with Br/Ag(100) in Section El Our conclusions are summarized 
in Section El 
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2 Self-consistent electrostatic adlayer model 



2.1 Lattice-gas model 

The adsorption of CI ions occurs at the fourfold hollow sites of the Ag(100) 
surface j2] > which form a square lattice as shown in Fig. ^ To approximate 
the equilibrium behavior of this system, we use a lattice-gas model, in which 
the lattice sites correspond to the adsorption sites. Mitchell et al. |T2] used an 
off-lattice model for the Br/Ag(100) system to show that the Br adsorbates 
spend most of the time near the four-fold hollow sites of the Ag(100) surface, 
thus justifying the lattice-gas treatment of halide adsorption. To describe 
the energy associated with a configuration of adsorbates on the surface, a 
grand-canonical effective Hamiltonian (21 El ECU El is used, 



where J2i<j is a sum over all pairs of sites, 0^- are the lateral interaction ener- 
gies between particles on the ith and jth lattice sites, measured in meV/pair, 
~p is the electrochemical potential, measured in meV/particle, and N = L 2 is 
the total number of lattice sites. The local occupation variable q is 1 if site 
i is occupied and otherwise. 

The long-range interactions, 0^, depend on the distance, r^, between 
ions i and j (measured in Ag(100) lattice spacing units, a = 2.889 A 0) as 



where the infinite value for = 1 indicates nearest-neighbor exclusion, and 
negative values of faj denote long-range repulsion. The coverage isotherms 
were simulated using a square LxL lattice with periodic boundary conditions 
to reduce finite-size effects. 

The electrochemical potential ~p is related to the bulk ionic concentration 
C and the electrode potential E (measured in mV). In the dilute-solution 
approximation, the relationship is 



H = - 4>ijCiCj 

i<3 



vJ2 



i=l 



(1) 






(3) 
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where /Z is an arbitrary constant, Co is a reference concentration (here taken 
to be 1 mM), and e is the elementary charge unit ^S]- The reference potential 
Eq is chosen sufficiently negative such that the coverage vanishes at Eq for all 
values of C used, and /Z has the sign convention that /Z > favors adsorption. 
The relationship between /J, C, and .E is discussed further in the Appendix. 

2.2 Lateral interaction energies 

When CI ions adsorb on the surface, a fraction of their charge is transfered 
through the external circuit. This fraction, je, is negative and is directly 
related to the average resident charge per ion, q = — (1 + j)e [TK] . This 
relationship is an approximation and is more valid as the potential at the 
adsorbate approaches the value of the potential in the solution. For the 
current system's ionic strength, this condition is only approximately satisfied 
and may be considered as a source of error. 

We have previously shown [3| that for Cl/Ag(100), the electrosorption 
valency 7 depends on the coverage 9, which is defined as 



In order to investigate such a dependence more thoroughly, we here propose 
a model with a coverage-dependent next-nearest-neighbor lateral interaction 
energy nnn , as well. This is motivated by two assumptions: that 7 is coverage 
dependent and that the major contribution to </> nnn is due to electrostatic 
dipole-dipole interactions. A simple electrostatic picture of the adlayer, in 
which an adsorbate's resident charge and its image charge form a dipole, 
suggests a relationship between the electrosorption valency and the dipole 
moment. If 7 is coverage dependent (7 = j(9)), then the resident charge is 
also coverage dependent (q = q(0)), and hence nnn , which is proportional to 
q 2 , is coverage dependent as well (0 nn n = 0nnn(#))- 

Assuming for simplicity that 7 depends linearly on 9, 



N 



9 = N- 1 Y,c l . 



(4) 



7 = 7o + 7i# > 



(5) 



the resident charge q becomes 




(6) 
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and nnn takes the form: 



nnn 



A{1 + 70 + liOf , 



(7) 



where A is a prefactor proportional to the square of a "dipole distance" . This 
effect is known as depolarization, and may under extreme conditions lead to 
the possibilty of a first-order phase transition ^7J ITHj 1 . If 7 is coverage 
independent (i.e., 71 =0), then nnn becomes coverage independent, as well. 

2.3 Monte Carlo method 

Equilibrium Monte Carlo (MC) simulations were used to measure the equi- 
librium coverage as a function of /Z, which was then converted to E using 
Eq. (jSJ). At each MC step a lattice site, i, was chosen randomly, and a change 
in its occupation variable Cj was attempted with a Metropolis acceptance 



where AH is the energy difference between the initial state and the proposed 
state. The value of nnn was updated in the simulations at each MC step, 
corresponding to the new proposed coverage. 

The energy difference AH was calculated from Eq. (Q) using two methods, 
a truncated sum of the contributions of neighboring occupied sites up to five 
lattice spacings away pD El E], and a mean- field-enhanced truncated sum 
up to three lattice spacings In the latter method, energy contributions 
from adparticles more than three lattice spacings away are calculated using 
a mean-field estimate as detailed in Ref. 

3 Experimental 

A detailed description of the experimental procedure is given in Ref. [HUE], 
and only a brief summary follows here. The Ag(100) single-crystal electrodes 
were chemically etched in cyanide solution, rinsed in Milli-Q water, and care- 
fully annealed in a hydrogen flame. They were then quickly transfered into 

1 For our model of CI adsorption on Ag(100) with the fitting parameters obtained in 
this paper, the first-order phase transition is not observed in the simultions, even for 
temperatures as unrealistically low as 1 K. The model does produce a first-order transition 
at room temperature for some parameters that do not fit the experimental data. 



probability [TU] : 



1Z = min 1 
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the electrochemical cell after cooling in a stream of Argon. Using the hang- 
ing meniscus technique, a platinum wire counter electrode, and a saturated 
calomel electrode as reference electrode, the electrochemical measurements 
were carried out at a temperature of (20 ± 1)°C. 

For the chronocoulometric experiments, the potential was set at initial 
values between —1.375 and —0.300 V vs SCE until adsorption equilibrium 
was established, and then stepped to the final potential of —1.400 V, where 
the chloride is completely desorbed from the surface. 

4 Results 

The equilibrium MC simulations were performed at a temperature of 17°C 
for various values of A, 70, and 71. The value of A was varied from to 
— 100 meV in steps of 5 meV, while 70 and 71 were varied from to —0.9 
in steps of 0.1, resulting in 2000 different simulations. The resulting sim- 
ulated isotherms were fitted to experimental data. This large number of 
simulations necessitated the choice of the relatively small system size with 
L = 32. The system displays a second-order phase transition from a disor- 
dered, low-density phase to an ordered c(2 x 2) phase [2] at an intermediate 
value of the electrochemical potential. For the steepest part of the isotherm 
(Fig. |2I), where critical slowing down due to the phase transition is impor- 
tant, the system was allowed to relax for a longer time to reach equilibrium 
than for the points further away from the phase transition. We checked for 
finite-size effects, and no significant difference was observed when comparing 
the resulting isotherms with isotherms simulated using larger system sizes 
of L = 64, 128, and 256. A small value of L can be used because 9 is not 
the order parameter corresponding to the c(2 x 2) phase, and so its fluctu- 
ations, which are proportional to dO/dp, only diverge logarithmically with 
L |2J 1201 121] ■ In addition, better statistics were collected close to the phase 
transition, using longer runs and more sampling than for points far from the 
phase transition, where the fluctuations are smaller. 

Using Eq. (j3J), the simulated isotherms were converted from the ~p scale 
to the E scale and then fit to the experimental isotherms. The fitting was 
done by varying the value of ~p to minimize \ 2 for each set of values for A, 
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7o, and 71. Here, 

^ ^ (^cxp ^sim) 
~2 _ cxp points / , 

degrees of freedom 

is the least-squares sum per degree of freedom, and ^exp and 9 S i m are the ex- 
perimental and simulated coverages, respectively, corresponding to the same 
value of E. Linear interpolation was used between the simulated data points 
to calculate 6 S i m for all values of E. 

From fits to 10 mM and 20 mM experiments, a grid of of values of % 2 , 
was collected. Three different models were studied, (i) A constant 7 (i.e., 
also constant n nn) model for which we only considered the simulations with 
71 = 0. (ii) A model in which only the line with 70 = —0.3, 71 = —0.3 
was considered. These values for 70 and 71 were estimated from values of 7 
obtained from the experimental data in Ref . [3] using the method in Ref . [22] ■ 
(hi) The third model studied was the coverage dependent model with nnn = 
^nnn(7(^)) ; where no constraints on the values of any of the parameters were 
imposed, and the whole parameter space was searched for a global minimum. 

A three-dimensional plot of the parameter space is shown in Fig. E^a), 
where the diameters of the symbols are proportional to 1/x 2 for fits to 
20 mM (circles) and 10 mM (squares) of the mean-field-enhanced simulations. 
The parameter space for the non-mean-field-enhanced method (not shown) is 
similar. Figure ED^a) shows the existence of several local minima and a good 
overlap between the minima for both concentrations. The position of the 
accepted global minimum is indicated by the arrows. Figure Efb) is a pro- 
jection onto the 70, 71 plane, which shows that the minima are concentrated 
within one region close to the 71 = plane, suggesting a relatively weak de- 
pendence of 7 on the coverage. Moreover, the x 2 values in the 71 = plane 
(model (i)) and for 70 = 7i = —0.3 (model (ii)) are significantly larger than 
the x 2 values for the accepted global minimum (model (hi)). Figures Efc) 
andEfd) show that while the magnitude of 70 decreases monotonically with 
increasing A, there is no general trend for 71 as a function of A. 

Due to the existence of several shallow minima of x 2 , and due to the 
limited resolution of the grid in parameter space, we list in Table 1 the best- 
fit parameter values for each model, along with all fits that have x 2 within 
10% of the best-fit value. One can see in Table 1 that there are two possible 
sets of parameters that fit to the 10 mM data for model (iii) with mean- 
field-enhanced simulations. To discriminate between these possible fits to the 
10 mM data we check if they are also possible fits to the 20 mM experimental 



7 



data. The first set (A = —90 meV) is the overall best fit for 10 mM. However, 
it fits only slightly better than the second set (A = —55 meV), and it is not 
a possible fit to the 20 mM experimental data. The only parameter set 
(A = —55 meV, 70 = —0.4, 71 = —0.2) that is a possible fit to both the 
10 mM and 20 mM experimental data, is the accepted best fit for model (iii). 
Using the same approach for the non-mean-field-enhanced method, models 
(i) and (ii), and concentration, results in a unique set of parameter values 
for each model, simulation method, and concentration. These accepted fits 
are summarized in Table 2. Notice in Table 2 that for the accepted global 
minimum of model (iii), and for both mean- field-enhanced and non- mean- 
field-enhanced methods, the values of 70 and 71 are reasonably close to the 
ones calculated from the experimental data j3] using the method of Ref . [22] • 
Also note that the fits for 10 mM are consistently better than the fits to the 
20 mM experiments. This might be due to violation of the dilute-solution 
limit at the higher concentration. This may also be the reason that the 
values of ~p , which are expected to be the same for both concentrations 
when fit to a single simulation, differ consistently by about 12 ± 7 meV 
between the concentrations. The values of A are consistently less negative 
for the longer-ranged lateral interactions of the mean-field-enhanced than for 
the non-mean-field method. This is not surprising. 

The quality of the fits is better for model (iii), as can be seen from the 
plots corresponding to each of the models (Figs.|21IH andEJ)- Models (ii) and 
(i) fit worse at either the lower- coverage part of the isotherm (model (ii)), 
see Fig. |U or at the upper part of the isotherm (model (i)), see Fig. El In 
contrast, the best-fit simulations in Fig. |2] fit the experiments well over the 
whole range of coverages. The plots for the non-mean-field-enhanced fits (not 
shown) are similar. A plot of 7 vs E is shown is Fig. Efa) for the best-fit 
values of 70 and 71 corresponding to the 10 mM experimental data for each 
of the three models considered. Also shown are the values of 7 obtained in 
Ref. jS| by the method of Ref. [221 • While model (ii) is expected to fit the 
Ref. [HI 122] values of 7, model (iii) also predicts a mean value of 7 close to, 
but slightly more negative than, the mean of the values of 7 from Ref. [3J 122] • 
The nnn values corresponding to the best fits of the three models are shown 
in Fig. Efb). 

Since the fits are sometimes better for the mean-field-enhanced method, 
and sometimes worse, the mean of the sets of parameters obtained by includ- 
ing and excluding the mean-field enhancement are reported here as our final 
results: A = —60 meV, 70 = —0.4, 7! = —0.2, and 7J = —325 meV. In the 



8 



far-field approximation used here, the lateral interaction energy between two 
parallel dipoles separated by a distance ay/2 is given by 



^=i^(i=^ 1+ * + ^ 2 =^(!) 2 ' (10) 

where p = \q\d is the dipole moment, d is an effective "dipolar distance," and 
the second and third equalities result from Eqs. (JTj) and ©, respectively. 
Consequently, 




For the A = —60 meV final result, the dipole moment ranges between 
p = 0.32 eA for 6 = and p = 0.26 eA for 9 = 0.5, and the "dipolar 
distance" d ~ 0.53 A. Density Functional Theory calculations of a model 
without water j2H| suggest a dipole moment of 0.46 eA, which is only about 
40% higher than the fitted value at low coverage. From naive physical intu- 
ition, the dipolar distance might be expected to be of the order of the ionic 
diameter of CI, about 3A. The value of d obtained here involves several as- 
sumptions. One main assumption used is that the dipole is a classical point 
dipole, which should be quite reasonable as the distance between neighbor- 
ing dipoles is about an order of magnitude larger than the obtained value 
for d. Moreover, we have assumed that the interaction is a classical electro- 
static dipole interaction, while there could be quantum effects at the shortest 
length scales. On the other hand, when viewed as a charge distribution, the 
dipole is expected to have a much smaller "dipolar distance" than naively 
expected |24j . Using the relationship between the dipole moment and the 
Helmholtz capacitance Ch [TE] , 

p=|?|d=^(l + 7), (12) 

and our values for the dipole moment and 7, we obtain a Ch that is in the 
range of 26 to 32 /iF/cm 2 . This range is comparable to the reported value of 
Ch for Ag(110) in chloride ion solution, Ch ~ 100 /zF/cm 2 [21], which would 
yield d = e /Cn ~ 0.1 A. This value of d is of the same order of magnitude 
as our value. 

Finally, although models (i) and (ii) are still possible, a self-consistent 
(</>nnn(7(#))) entirely electrostatic model not only fits better to the experi- 
mental data but also, with no constraints imposed, predicts values of 7 that 
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are reasonably compatible with those obtained from the experimental data 
using the method of Ref. [22] . 

5 Comparison to Br/Ag(100) 

Br electrosorption on the same substrate Ag(100) displays different char- 
acteristics. The electrosorption valency and the next-nearest-neighbor lat- 
eral interaction energy are not only more negative, 7 = —0.71 ± 0.01 and 
0nnn = —21 ±2 meV, but also independent of the coverage [3]. In other words, 
while mutual depolarization JH] is present for Cl/Ag(100), this effect is not 
significant for Br/Ag(100). Simulations of Br/Ag(100) using 7 = 7o + 7i#, 
when fit to three experimental data sets, yielded 71 ~ |3 a . Moreover, the 
value obtained for 7 (—0.71), is consistent with previous results from both 
simulations and experimental analysis |2|. Since CI is more electronegative 
than Br, it is expected to have a less negative 7 or a more negative resident 
charge, consistent with our results. 

6 Conclusion 

A discrepancy between the values of the electrosorption valency 7 obtained 
from fitting simulations to experiments, and the value obtained from the 
analysis [22] of experimental data was reported in a previous study by Abou 
Hamad et al. [3J. This discrepancy suggested that the long-range interactions 
are dominated by electrostatic dipole/dipole effects for CI electrosorption on 
the Ag(100) single crystal surface. 

In this work, a large set of MC simulations (2000 simulations) over a grid 
in parameter space was fitted to two sets of chronocoulometry experimental 
data for different concentrations. The existence of local minima of x 2 in pa- 
rameter space suggests alternative models. To within the resolution of our 
grid and the accuracy of the experimental data, we have shown that while 
other, simpler, models are still possible, a purely electrostatic model can be 
used to describe the Cl/Ag(100) system. It also predicts an electrosorption 
valency that is compatible with the value obtained through direct experimen- 
tal data analysis by the method of Ref. [221- Additional sets of experimental 
data, with different concentrations around 10 and 20 mM of the adsorbate 
ion, along with a finer grid in parameter space would enable a more decisive 
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determination of the most appropriate model and its parameter values. 
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In this appendix we discuss the relationship between the electrochemical 
potential /Z and the electrode potential E for a coverage or field dependent 
electrosorption valency 7. 

We define the electrosorption valency as [To] 



where <tm is the charge on the metal, 6 is the coverage, and /1 is the chemical 
potential of the solution (here taken as the dilute-solution approximation, 
\i = k^T\ia(-^)). Considering an auxiliary thermodynamic potential X re- 
lated to the surface tension a as 



Appendix 




(13) 



X = a + Ojj, 



(14) 



and using the Lippman equation, 



dX = idd6 - a M dE , 



(15) 



we have 




(16) 



Consequently, one obtains the Maxwell relation, [16 




(17) 
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The relation, 7 = (f§J , gives the electrosorption valency as the change 
in \i necessary to keep constant under a change in E. But to keep 6 
constant, the electrochemical potential ~p must remain constant, so that 

It is easy to show that Eq. (JHJ) satisfies this requirement for 7. If 

Je 



then 



and 



/I = 7Z + fc B Tln (— ) - e ^ 7 (S')dS / , 

pi = -p+e / i(E')&E' + const. , (19) 



dfi\ ( d[i 
dE) n ~ [dE 



= 1(E) ■ (20) 
Generalization to multiple adsorbates is straightforward. 
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Figure 1: CI (larger, dark gray, spheres) adsorbed at the 4-fold hollow sites of 
the (100) surface of Ag (smaller, lighter gray, spheres). The grid frame corre- 
sponds to the lattice of adsorption sites. The figure is drawn approximately 
to scale. 
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E/VvsSCE 

Figure 2: Best fit of mean-field-enhanced simulation to experimental 10 mM 
and 20 mM coverage isotherms: A = —55 meV, 70 = —0.4, and 71 = —0.2 
(model (iii)). L = 32. 
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Figure 3: A plot of 1/x 2 in the three-dimensional parameter space, obtained from fits of mean-field-enhanced 
simulations to the 10 mM experimental data (squares) and the 20 mM experimental data (circles). The 
diameters of the symbols are proportional to 1/x 2 - (a) is a three dimensional view, (b) is a projection onto 
the 70, 71 plane, and (c) and (d) are projections onto the 7 , A plane and the 71, A, plane respectively. The 
parameter space for the non-mean-field-enhanced method (not shown) is similar. 




-0.25 



E/VvsSCE 

Figure 4: Simulated constrained best fit with the values of 70 = —0.3 and 
71 = —0.3 (model (ii)) for the 10 and 20 mM experimental data. Mean-field- 
enhanced simulations are shown here. L = 32. 
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E/VvsSCE 

Figure 5: Simulated constrained best fit with a constant 7 (model (i)) for 
the 10 and 20 mM experimental data. Mean-field-enhanced simulations are 
shown here. L = 32. 
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q-O Y = -0.4, y = 0, model (i) 
T Ref. [3,19] y 



i,25 1 1 1 1 1 1 1 1 

-1.25 -1 -0.75 -0.5 -0.25 

E / V vs SCE 




E / V vs SCE 

Figure 6: (a) The electrosorption valency 7 vs E for the values of 70 and 
71 obtained from the best fits for models (i), (ii), and (iii) to the 10 mM 
experimental data. The values of 7 obtained in Ref. [3] by the method of 
Ref. |22| are also shown, (b) Corresponding plots of nnn vs E for the three 
models. 
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Table 1: Fitting parameters for CI adsorption on a Ag(100) single-crystal 
surface, for models (i), (ii), and (iii), with and without mean-field interac- 
tions. Possible fits to within 10% of the best x 2 f° r each model are also 
shown. A and ~p Q are in units of meV. 



Model 


Mean-field enhanced 


10 mM 


20 mM 


(i) Const. 7 


A 


7o 


7i 


^0 


X 2 x 10 5 


A 


7o 


7i 


^0 


X 2 x 10 5 


-25 


-0.4 





-330 


5.664 


-60 


-0.5 





-390 


11.927 


-5 


-0.3 





-281 


5.769 


-55 


-0.5 





-396 


11.942 












-25 


-0.4 





-340 


12.079 












-20 


-0.4 





-349 


12.514 


(ii) Ref. 3, 221 7 


-30 


-0.3 


-0.3 


-266 


6.068 


-25 


-0.3 


-0.3 


-285 


8.037 


(iii) 0„ nn (7W) 


-90 


-0.5 


-0.1 


-371 


1.815 


-55 


-0.4 


-0.2 


-328 


3.340 


-55 


-0.4 


-0.2 


-318 


1.976 












(i) Const. 7 


No mean-field enhancement 


-30 


-0.4 





-330 


5.369 


-70 


-0.5 





-393 


11.246 












-30 


-0.4 





-341 


11.330 












-65 


-0.5 





-398 


11.350 












-25 


-0.4 





-348 


11.606 


(ii) Ref. 3, 221 7 


-35 


-0.3 


-0.3 


-268 


6.086 


-30 


-0.3 


-0.3 


-286 


7.462 


(iii) 0nnn(7(#)) 


-45 


-0.4 


-0.1 


-325 


2.310 


-65 


-0.4 


-0.2 


-330 


3.112 


-65 


-0.4 


-0.2 


-320 


2.311 













Table 2: The best-fit parameters for CI adsorption on a Ag(100) single- 
crystal surface, for models (i), (ii), and (iii), with and without mean-field 
interactions. A and ~p are in units of meV. 



Model 


Mean-field enhanced 


10 mM 


20 mM 


(i) Const. 7 


A 


7o 


7i 




X 2 x 10 s 


A 


7o 


7i 


Mo 


X 2 x 10 b 


-25 


-0.4 





-330 


5.664 


-25 


-0.4 





-340 


12.079 


(ii) Ref. 3, 22] 7 


-30 


-0.3 


-0.3 


-266 


6.068 


-25 


-0.3 


-0.3 


-285 


8.037 


(iii) nrm (7W) 


-55 


-0.4 


-0.2 


-318 


1.976 


-55 


-0.4 


-0.2 


-328 


3.340 


(i) Const. 7 


No mean-field enhancement 


-30 


-0.4 





-330 


5.369 


-30 


-0.4 





-341 


11.330 


(ii) Ref. 3,1221 7 


-35 


-0.3 


-0.3 


-268 


6.086 


-30 


-0.3 


-0.3 


-286 


7.462 


(iii) <?Wi(7(0)) 


-65 


-0.4 


-0.2 


-325 


2.310 


-65 


-0.4 


-0.2 


-330 


3.112 
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